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Рассмотрены математические методы оценки корреляции неисправностей локомотивов. Представлен способ оптимизации мно- 
жества диагностируемых параметров. Показаны возможности прогнозирования технического состояния диагностируемого 
объекта. 


Существующие на сегодня ш ний день методы и 
средства технического диагностирования локомо- 
тивов позволяют оценить их техническое состояние 
после возвращения в депо из поездки. Однако зна- 
чительный объём информации, получаемой в про- 
цессе движения, не используется. Источником по- 
добной информации может служить разработанная 
бортовая система диагностирования. При опти- 
мальном использовании дополнительной инфор- 
мации становится возможным повышение досто- 
верности оценки технического состояния электро- 
возов. Этого можно достичь путём увеличения ко- 
личества обрабатываемых данных, установления за- 
висимостей в причинах изменения технического 
состояния и на основе полученной информации 
прогнозировать возможные неисправности. На се- 
годняшний день актуальность поставленной задачи 
возрастает в связи с общим старением локомотив- 
ного парка, а также необходимостью повышения 
безопасности эксплуатации подвижного состава. 

Для контроля работоспособности электропо- 
движного состава требуется подтверждение нор- 
мального функционирования всех его элементов. 
Решение этой проблемы усложняется необходимо- 
стью проведения значительных объёмов и количе- 
ства контрольных операций [1]. Использование 
информации бортовой системы диагностирования, 
полученной в ходе проведения диагностических 
тестов электровоза, позволит перейти от плановых 
ремонтов к ремонтам по фактическому техниче- 
скому состоянию, снизить затраты на техническое 
обслуживание и профилактику. 

Решение поставленной задачи требует опреде- 
ления причин изменения технического состояния. 
Обычно для определения степени связанности двух 
событий используют значение коэффициента кор- 
реляции. Однако проблема, указывает ли получен- 
ная оценка коэффициента корреляции на зависи- 
мость двух событий с требуемой достоверностью 
или нет, остаётся и требует решения. 

Сформулируем задачу следующим образом: 

У (к) = КХ(сІ), 

где У- множество оценок технического состояния 
объектов контроля к, X- множество диагностируе- 
мых параметров сі. 


Требуется найти оператор отображения К мно- 
жества Хна множество У. Ставится задача найти та- 
кой оператор К, при котором скалярное произведе- 
ние множества Ун У*, должно стремиться к 1 . Здесь 
У* действительное значение оценок технического 
состояния. Элементы оператора К являются коэф- 
фициентами корреляции между элементами мно- 
жеств Хи У. 

Корреляционный момент двух случайных вели- 
чин х, у определяется по следующей формуле: 

со оо 

К ху = | (х-ш г Щх)йС | (у -т у )/ 2 (у)с1у. 

—00 —00 

где т„ т у - моменты выборок хи у, соответствен- 
но;/^),^^) - функции плотности распределения 
вероятности х и у, соответственно. 

Пусть (Хі , у ! ) , (х 2 ,^ 2 ) , . . . , (х„ ,у я ) - выборка из и на- 
блюдений пары случайных переменных ( X , У), ко- 
торые имеют совместное двумерное нормальное 
распределение с коэффициентом корреляции д. 
Обозначим х=Ух/п, Ѵ=У.у/п. Коэффициент корре- 
ляции выборки г определяется как 

г ху = а/^/іВС) , ( 1 ) 

где 

А = X '(х ~ х )(Уі - у ')} = -пху, 

1 1 

В =^(х-х) 2 =^х 2 -пх 2 , 

1 1 

С = Х(Х -у ) 2 =Т^Уі~пу 2 - 
1 1 

Выборочное распределение коэффициента 
корреляции выборки имеет плотность распределе- 
ния вероятности в точке г, заданную в виде 

/„ (г, д) = (1 - д 2 Г ,П (\ ~ г 2 У~* )І2 а(п, г, д), 

-\<г<\, 

где а(п, г, д) = ^ Г 2 {(и + 5 - 1)/2} (2дг) 1 /я !, 

Г{} — гамма-функция, п - число наблюдений, 
д - коэффициент корреляции двух переменных 
(Х,У), г - коэффициент корреляции выборки ( 1 ). 


155 




Известия Томского политехнического университета. 2006. Т. 309. № 1 


Для характеристики связи между величинами 
( X , У) в чистом виде перейдём от момента К„ к безраз- 
мерной характеристике коэффициента корреляции 


где <у„ сг ѵ . - средние квадратичные отклонения ве- 
личин X, У. 

Теперь следует определить оптимальное множе- 
ство диагностируемых параметров, на котором до- 
стоверность оценки технического состояния до- 
стигает максимума. Оптимизация этого множества 
будет производиться путём исключения из него не- 
связанных данных, критерием отбора будет слу- 
жить некоторое условие связанности данных. 

Выдвинем гипотезу Н о том, что полученный 
коэффициент корреляции действительно позволя- 
ет судить о связи имеющихся данных. Когда 
необходимо проверить, указывает ли наблюдаемое 
значение г на действительную коррелированность 
данных, то соответствующая нулевая гипотеза Н 
имеет вид: 

Н: <7=0. 

Достаточно большое по абсолютной величине г 
будет стремиться опровергнуть нулевую гипотезу [2]. 

В рассматриваемой задаче обнаружения зависи- 
мости между различными типами неисправностей 
и измерений опровержение этой гипотезы приве- 
дет к подтверждению существования указанной за- 
висимости. На этот вопрос можно ответить с помо- 
щью преобразования 

' = П=Т>- 2)- (2) 

у/1 г *у 

Итак, если корреляция двух каких-либо собы- 
тий превышает допустимый порог, то полагаем их 
взаимосвязанными, т.е., например, возникновение 
неисправности А влечет за собой появление неис- 
правности В, или же обе эти неисправности явля- 
ются следствием, ранее произошедшего события С. 
Это приведёт в перспективе к более точному прог- 
нозированию событий и выявлению связей между 
возникновением неисправностей в одном агрегате 
и сбоев в работе другого. Что в свою очередь увели- 
чит эффективность диагностирования. 

При нулевой гипотезе выборочное распределе- 
ние статистики неисправностей есть распределение 
Стьюдента с я-2 степенями свободы. Большие абсо- 
лютные значения г отвечают большим абсолютным 
значениям /, а поскольку выборочное распределение 
1 симметрично относительно точки 0, то становится 
возможным определить уровень значимости 8Ь: 

8Ь=Р(Т>\і\)+Р(Т<-Щ)=2Р(Щ), 
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где Т подчиняется распределению Стьюдента с я-2 
степенями свободы, а 1 выбирается в соответствии с 
(2) по выборочному коэффициенту корреляции г [2]. 

Таким образом, был получен уровень зависимо- 
сти между событиями, который позволит получить 
однозначный ответ на вопрос, что может служить 
причиной отказа. В контексте задачи диагностиро- 
вания ими могут быть: неисправности, включая 
время и место их возникновения, наличие внеш- 
них воздействий и т.д. 

Для решения задачи оптимизации множества 
диагностируемых параметров следует исключить 
все данные, значения коэффициентов корреляции 
которых ниже уровня значимости. Это позволит 
избавиться от лишнего анализа несвязанных собы- 
тий. Оптимизация представляет собой итерацион- 
ный процесс, на каждом шаге которого происходит 
сравнение и оценивание параметров. 

Пусть, например, X- множество отказов и из- 
мерений, а У множество отказов, т.е. У- подмно- 
жество X. Найдём такие X * - исключаемые элемен- 
ты множества X, 

х*еХ*, і = 1, 2Ѵ, 

для которых 

Исключив их из рассматриваемого множества, по- 
лучим новые оптимизированные множества Хц и У 0 . 

В качестве примера возьмём выборочный коэф- 
фициент кор реляции между двумя оцениваемыми 
параметрами г= — 0,63. Соответствующее значение 1, 
вычисленное по формуле (2) для набора из 20 наблю- 
дений равно — 3,43. Уровень значимости составляет 
2Р(Г 20 >5,44). Значению 1=-3,43 в соответствии с ра- 
спределением Стьюдента отвечает значение 8Ь, ме- 
ньшее 0,01, что значимо. Нулевая гипотеза отклоне- 
на, а существование корреляции установлено. 

Применение соответствующих математических 
методов даёт возможность прогнозировать техни- 
ческое состояние диагностируемого объекта и по- 
зволяет сократить множество элементов объекта 
диагностирования и, как следствие, повысить точ- 
ность в определении неисправностей, что приведёт 
к повышению безопасности движения. 

Заключение 

На практике учёт коэффициента корреляции 
неисправностей диагностируемого оборудования 
позволит повысить достоверность диагностиче- 
ской информации, улучшить качество диагности- 
рования агрегатов и узлов, выявить связи между 
различными событиями. Это позволит продлить 
срок работы локомотива между ремонтами, сни- 
зить степень износа оборудования, экономить ма- 
териальные и трудовые ресурсы. 
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